Yam genomics supports West Africa as a major cradle of crop domestication.

While there has been progress in our understanding of the origin and history of agriculture in sub-Saharan Africa, a unified perspective is still lacking on where and how major crops were domesticated in the region. Here, we investigated the domestication of African yam (Dioscorea rotundata), a key crop in early African agriculture. Using whole-genome resequencing and statistical models, we show that cultivated yam was domesticated from a forest species. We infer that the expansion of African yam agriculture started in the Niger River basin. This result, alongside with the origins of African rice and pearl millet, supports the hypothesis that the vicinity of the Niger River was a major cradle of African agriculture.

Assessing the Response of Small RNA Populations to Allopolyploidy Using Resynthesized Brassica napus Allotetraploids.

Allopolyploidy, combining interspecific hybridization with whole genome duplication, has had significant impact on plant evolution. Its evolutionary success is related to the rapid and profound genome reorganizations that allow neoallopolyploids to form and adapt. Nevertheless, how neoallopolyploid genomes adapt to regulate their expression remains poorly understood. The hypothesis of a major role for small noncoding RNAs (sRNAs) in mediating the transcriptional response of neoallopolyploid genomes has progressively emerged. Generally, 21-nt sRNAs mediate posttranscriptional gene silencing by mRNA cleavage, whereas 24-nt sRNAs repress transcription (transcriptional gene silencing) through epigenetic modifications. Here, we characterize the global response of sRNAs to allopolyploidy in Brassica, using three independently resynthesized Brassica napus allotetraploids originating from crosses between diploid Brassica oleracea and Brassica rapa accessions, surveyed at two different generations in comparison with their diploid progenitors. Our results suggest an immediate but transient response of specific sRNA populations to allopolyploidy. These sRNA populations mainly target noncoding components of the genome but also target the transcriptional regulation of genes involved in response to stresses and in metabolism; this suggests a broad role in adapting to allopolyploidy. We finally identify the early accumulation of both 21- and 24-nt sRNAs involved in regulating the same targets, supporting a posttranscriptional gene silencing to transcriptional gene silencing shift at the first stages of the neoallopolyploid formation. We propose that reorganization of sRNA production is an early response to allopolyploidy in order to control the transcriptional reactivation of various noncoding elements and stress-related genes, thus ensuring genome stability during the first steps of neoallopolyploid formation.

Systematic comparison of small RNA library preparation protocols for next-generation sequencing.

BACKGROUND: Next-generation sequencing technologies have revolutionized the study of small RNAs (sRNAs) on a genome-wide scale. However, classical sRNA library preparation methods introduce serious bias, mainly during adapter ligation steps. Several types of sRNA including plant microRNAs (miRNA), piwi-interacting RNAs (piRNA) in insects, nematodes and mammals, and small interfering RNAs (siRNA) in insects and plants contain a 2'-O-methyl (2'-OMe) modification at their 3' terminal nucleotide. This inhibits 3' adapter ligation and makes library preparation particularly challenging. To reduce bias, the NEBNext kit (New England Biolabs) uses polyethylene glycol (PEG), the NEXTflex V2 kit (BIOO Scientific) uses both randomised adapters and PEG, and the novel SMARTer (Clontech) and CATS (Diagenode) kits avoid ligation altogether. Here we compared these methods with Illumina's classical TruSeq protocol regarding the detection of normal and 2' OMe RNAs. In addition, we modified the TruSeq and NEXTflex protocols to identify conditions that improve performance. RESULTS: Among the five kits tested with their respective standard protocols, the SMARTer and CATS kits had the lowest levels of bias but also had a strong formation of side products, and as a result performed relatively poorly with biological samples; NEXTflex detected the largest numbers of different miRNAs. The use of a novel type of randomised adapters called MidRand-Like (MRL) adapters and PEG improved the detection of 2' OMe RNAs both in the TruSeq as well as in the NEXTflex protocol. CONCLUSIONS: While it is commonly accepted that biases in sRNA library preparation protocols are mainly due to adapter ligation steps, the ligation-free protocols were not the best performing methods. Our modified versions of the TruSeq and NEXTflex protocols provide an improved tool for the study of 2' OMe RNAs.

Molecular basis of African yam domestication: analyses of selection point to root development, starch biosynthesis, and photosynthesis related genes.

BACKGROUND: After cereals, root and tuber crops are the main source of starch in the human diet. Starch biosynthesis was certainly a significant target for selection during the domestication of these crops. But domestication of these root and tubers crops is also associated with gigantism of storage organs and changes of habitat. RESULTS: We studied here, the molecular basis of domestication in African yam, Dioscorea rotundata. The genomic diversity in the cultivated species is roughly 30% less important than its wild relatives. Two percent of all the genes studied showed evidences of selection. Two genes associated with the earliest stages of starch biosynthesis and storage, the sucrose synthase 4 and the sucrose-phosphate synthase 1 showed evidence of selection. An adventitious root development gene, a SCARECROW-LIKE gene was also selected during yam domestication. Significant selection for genes associated with photosynthesis and phototropism were associated with wild to cultivated change of habitat. If the wild species grow as vines in the shade of their tree tutors, cultivated yam grows in full light in open fields. CONCLUSIONS: Major rewiring of aerial development and adaptation for efficient photosynthesis in full light characterized yam domestication.

Polyploidy and interspecific hybridization: partners for adaptation, speciation and evolution in plants.

Background: Polyploidy or whole-genome duplication is now recognized as being present in almost all lineages of higher plants, with multiple rounds of polyploidy occurring in most extant species. The ancient evolutionary events have been identified through genome sequence analysis, while recent hybridization events are found in about half of the world's crops and wild species. Building from this new paradigm for understanding plant evolution, the papers in this Special Issue address questions about polyploidy in ecology, adaptation, reproduction and speciation of wild and cultivated plants from diverse ecosystems. Other papers, including this review, consider genomic aspects of polyploidy. Approaches: Discovery of the evolutionary consequences of new, evolutionarily recent and ancient polyploidy requires a range of approaches. Large-scale studies of both single species and whole ecosystems, with hundreds to tens of thousands of individuals, sometimes involving 'garden' or transplant experiments, are important for studying adaptation. Molecular studies of genomes are needed to measure diversity in genotypes, showing ancestors, the nature and number of polyploidy and backcross events that have occurred, and allowing analysis of gene expression and transposable element activation. Speciation events and the impact of reticulate evolution require comprehensive phylogenetic analyses and can be assisted by resynthesis of hybrids. In this Special Issue, we include studies ranging in scope from experimental and genomic, through ecological to more theoretical. Conclusions: The success of polyploidy, displacing the diploid ancestors of almost all plants, is well illustrated by the huge angiosperm diversity that is assumed to originate from recurrent polyploidization events. Strikingly, polyploidization often occurred prior to or simultaneously with major evolutionary transitions and adaptive radiation of species, supporting the concept that polyploidy plays a predominant role in bursts of adaptive speciation. Polyploidy results in immediate genetic redundancy and represents, with the emergence of new gene functions, an important source of novelty. Along with recombination, gene mutation, transposon activity and chromosomal rearrangement, polyploidy and whole-genome duplication act as drivers of evolution and divergence in plant behaviour and gene function, enabling diversification, speciation and hence plant evolution.

The bright side of transposons in crop evolution.

The past decades have revealed an unexpected yet prominent role of so-called 'junk DNA' in the regulation of gene expression, thereby challenging our view of the mechanisms underlying phenotypic evolution. In particular, several mechanisms through which transposable elements (TEs) participate in functional genome diversity have been depicted, bringing to light the 'TEs bright side'. However, the relative contribution of those mechanisms and, more generally, the importance of TE-based polymorphisms on past and present phenotypic variation in crops species remain poorly understood. Here, we review current knowledge on both issues, and discuss how analyses of massively parallel sequencing data combined with statistical methodologies and functional validations will help unravelling the impact of TEs on crop evolution in a near future.

Allopolyploidy has a moderate impact on restructuring at three contrasting transposable element insertion sites in resynthesized Brassica napus allotetraploids.

The role played by whole-genome duplication (WGD) in evolution and adaptation is particularly well illustrated in allopolyploids, where WGD is concomitant with interspecific hybridization. This 'Genome Shock', usually accompanied by structural and functional modifications, has been associated with the activation of transposable elements (TEs). However, the impact of allopolyploidy on TEs has been studied in only a few polyploid species, and not in Brassica, which has been marked by recurrent polyploidy events. Here, we developed sequence-specific amplification polymorphism (SSAP) markers for three contrasting TEs, and compared profiles between resynthesized Brassica napus allotetraploids and their diploid Brassica progenitors. To evaluate restructuring at TE insertion sites, we scored changes in SSAP profiles and analysed a large set of differentially amplified SSAP bands. No massive structural changes associated with the three TEs surveyed were detected. However, several transposition events, specific to the youngest TE originating from the B. oleracea genome, were identified. Our study supports the hypothesis that TE responses to allopolyploidy are highly specific. The changes observed in SSAP profiles lead us to hypothesize that they may partly result from changes in DNA methylation, questioning the role of epigenetics during the formation of a new allopolyploid genome.

BraSto, a Stowaway MITE from Brassica: recently active copies preferentially accumulate in the gene space.

We characterized a Brassica miniature inverted repeat transposable element (MITE) from the Stowaway superfamily, designated BraSto (Bra ssica Sto waway). BraSto copy number was assessed using real-time quantitative PCR in the two diploid species B. rapa (genome A) and B. oleracea (genome C) and the corresponding allotetraploid species B. napus (genome AC). Phylogenetic relationships among a set of 131 BraSto copies were then analyzed. BraSto appears to have been only moderately amplified in the Brassica genome and was still active recently with marks of proliferation in both diploid Brassica species, which diverged 3.75 million years ago, but also in the allotetraploid species after reuniting of the two differentiated genomes. We characterized insertion sites for low-divergence BraSto copies among the gene space of the B. rapa genome using bioinformatics approaches. For BraSto copies localized nearby or within genes, we observed frequent associations of BraSto with putative promoters and regulatory regions of genes, but exclusion from coding regions. In addition, BraSto was significantly similar to several Brassica expressed sequence tags (ESTs), including stress-induced ESTs. We also demonstrated the enrichment of BraSto sequences in binding sites for transcription factors and other regulatory elements. Our results lead to the question of a role for BraSto in the regulation of gene expression: this putative role, if further confirmed experimentally, would help to obtain a new insight into the significance of MITEs in the functional plant genome.

Analysis of gene expression in resynthesized Brassica napus allotetraploids: transcriptional changes do not explain differential protein regulation.

Polyploidy, or whole genome duplication, is a major evolutionary process that has shaped eukaryotic genomes, notably those of flowering plants. The mechanisms underlying the regulation of, and sharing of functions between, the duplicated genes originating from polyploidy events, which lead to novel phenotypes, remain to be elucidated. A previous comparative proteomic study identified 360 proteins that were differentially regulated between the diploid Brassica progenitors and their synthetic allotetraploid derivatives. For 102 of these proteins, using the same resynthesized Brassica napus allotetraploids, we assayed the accumulation of the transcripts of the corresponding genes. We compared transcript levels quantified in the synthetic allotetraploids with the mid-parent expression values. Although all of the genes surveyed encoded nonadditive proteins, we found that two-thirds of them had additive transcript levels, indicating that most of the differential protein regulation is not explained by transcriptional changes. Our data suggest that differential protein regulation is mainly governed by post-transcriptional modifications. Summarizing available data from transcriptomic studies of other synthetic allopolyploid models, we describe the general trends of transcript regulation in an allopolyploid genome and discuss putative underlying molecular mechanisms, with particular emphasis on the small RNA pathway for the post-transcriptional control of gene expression.

Impact of transposable elements on the organization and function of allopolyploid genomes.

Transposable elements (TEs) represent an important fraction of plant genomes and are likely to play a pivotal role in fuelling genome reorganization and functional changes following allopolyploidization. Various processes associated with allopolyploidy (i.e. genetic redundancy, bottlenecks during the formation of allopolyploids or genome shock following genome merging) may allow accumulation of TE insertions. Our objective in carrying out a survey of the literature and a comparative analysis across different allopolyploid systems is to shed light on the structural, epigenetic and functional modifications driven by TEs during allopolyploidization and subsequent diploidization. The available evidence indicates that TE proliferation in the short or the long term after allopolyploidization may be restricted to a few TEs, in specific polyploid systems. By contrast, data indicate major structural changes in the TE genome fraction immediately after allopolyploidization, mainly through losses of TE sequences as a result of recombination. Emerging evidence also suggests that TEs are targeted by substantial epigenetic changes, which may impact gene expression and genome stability. Furthermore, TEs may directly or indirectly support the evolution of new functionalities in allopolyploids during diploidization. All data stress allopolyploidization as a shock associated with drastic genome reorganization. Mechanisms controlling TEs during allopolyploidization as well as their impact on diploidization are discussed.

Contrasting evolutionary patterns and target specificities among three Tourist-like MITE families in the maize genome.

Miniature inverted-repeat transposable elements (MITEs) are short, non autonomous DNA elements that are widespread and abundant in plant genomes. The high sequence and size conservation observed in many MITE families suggest that they have spread recently throughout their respective host genomes. Here we present a maize genome wide analysis of three Tourist-like MITE families, mPIF, and two previously uncharacterized families, ZmV1 and Zead8. We undertook a bioinformatic analysis of MITE insertion sites, developed methyl-sensitive transposon display (M-STD) assays to estimate the associated level of CpG methylation at MITE flanking regions, and conducted a population genetics approach to investigate MITE patterns of expansion. Our results reveal that the three MITE families insert into genomic regions that present specific molecular features: they are preferentially AT rich, present low level of cytosine methylation as compared to the LTR retrotransposon Grande, and target site duplications are flanked by large and conserved palindromic sequences. Moreover, the analysis of MITE distances from predicted genes shows that 73% of 263 copies are inserted at less than 5 kb from the nearest predicted gene, and copies from Zead8 family are significantly more abundant upstream of genes. By employing a population genetic approach we identified contrasting patterns of expansion among the three MITE families. All elements seem to have inserted roughly 1 million years ago but ZmV1 and Zead8 families present evidences for activity of several master copies within the last 0.4 Mya.

The CACTA transposon Bot1 played a major role in Brassica genome divergence and gene proliferation.

We isolated and characterized a Brassica C genome-specific CACTA element, which was designated Bot1 (Brassica oleracea transposon 1). After analysing phylogenetic relationships, copy numbers and sequence similarity of Bot1 and Bot1 analogues in B. oleracea (C genome) versus Brassica rapa (A genome), we concluded that Bot1 has encountered several rounds of amplification in the oleracea genome only, and has played a major role in the recent rapa and oleracea genome divergence. We performed in silico analyses of the genomic organization and internal structure of Bot1, and established which segment of Bot1 is C-genome specific. Our work reports a fully characterized Brassica repetitive sequence that can distinguish the Brassica A and C chromosomes in the allotetraploid Brassica napus, by fluorescent in situ hybridization. We demonstrated that Bot1 carries a host S locus-associated SLL3 gene copy. We speculate that Bot1 was involved in the proliferation of SLL3 around the Brassica genome. The present study reinforces the assumption that transposons are a major driver of genome and gene evolution in higher plants.

Repetitive sequence-derived markers tag centromeres and telomeres and provide insights into chromosome evolution in Brassica napus.

Centromeres and telomeres are obvious markers on chromosomes but their location on genetic maps is difficult to determine, which hampers many basic and applied research programmes. In this study, we used the characteristic distribution of five Brassica repeated sequences to generate physically anchored molecular markers tentatively tagging Brassica centromeres (84 markers) and telomeres (31 markers). These markers were mapped to the existing oilseed rape genetic map. Clusters of centromere-related loci were observed on 14 linkage groups; in addition to previous reports, we could thus provide information about the most likely position of centromeres on 17 of the 19 B. napus linkage groups. The location of centromeres on linkage groups usually matches their position on chromosomes and coincides with sites of evolutionary breakage between chromosomes. Most telomere sequence-derived markers mapped interstitially or in the proximity of centromeres; this result echoes previous reports on many eukaryote genomes and may reflect different forms of chromosome evolution. Seven telomere sequence-derived markers were located at the outermost positions of seven linkage groups and therefore probably tagged telomeres.

Differential regulation of gene products in newly synthesized Brassica napus allotetraploids is not related to protein function nor subcellular localization.

BACKGROUND: Allopolyploidy is a preeminent process in plant evolution that results from the merger of distinct genomes in a common nucleus via inter-specific hybridization. Allopolyploid formation is usually related to genome-wide structural and functional changes though the underlying mechanisms operating during this "genomic shock" still remain poorly known. The aim of the present study was to investigate the modifications occurring at the proteomic level following an allopolyploidization event and to determine whether these changes are related to functional properties of the proteins. In a previous report, we applied comparative proteomics to synthetic amphiploids of Brassica napus and to its diploid progenitors B. rapa and B. oleracea. Although several hundred polypeptides displayed additivity (i.e. mid-parent values) in the amphiploids, many of them showed non-additivity. Here, we report the in silico functional characterization of the "non-additive" proteins (the ones with a non-additive pattern of regulation) in synthetic B. napus. RESULTS: The complete set of non-additive proteins (335 in the stem and 205 in the root), as well as a subset of additive polypeptides (200 per organ), was identified by mass spectrometry. Several protein isoforms were found, and most of them (approximately 55%) displayed "different" or "opposite" patterns of regulation in the amphiploids, i.e. isoforms of the same protein showing both up-regulation and down-regulation in the synthetic B. napus compared to the mid-parent value. Components of protein complexes were identified of which approximately 50% also displayed "different" or "opposite" patterns of regulation in the allotetraploids. In silico functional categorization of the identified proteins was carried out, and showed that neither functional category nor metabolic pathway were systematically affected by non-additivity in the synthetic amphiploids. In addition, no subcellular compartment was found to be over- or under-represented among the proteins displaying non-additive values in the allopolyploids. CONCLUSION: Protein identification showed that functionally related polypeptides (isoforms and complex subunits) could be differentially regulated in synthetic B. napus in comparison to its diploid progenitors while such proteins are usually expected to display co-regulation. The genetic redundancy within an allopolyploid could explain why functionally related proteins could display imbalanced levels of expression. No functional category, no metabolic pathway and no subcellular localization was found to be over- or under-represented within non-additive polypeptides, suggesting that the differential regulation of gene products was not related to functional properties of the proteins. Thus, at the protein level, there is no evidence for the "genomic shock" expected in neo-polyploids and the overall topology of protein networks and metabolic pathways is conserved in synthetic allotetraploids of B. napus in comparison to its diploid progenitors B. rapa and B. oleracea.

The genomic organization of retrotransposons in Brassica oleracea.

We have investigated the copy numbers and genomic organization of five representative reverse transcriptase domains from retrotransposons in Brassica oleracea. Two non-homologous Pseudoviridae (Ty1/copia-like) elements, two Metaviridae (Ty3/gypsy-like) elements (one related to the Athila family) and one Retroposinae (LINE) element were hybridized to a gridded BAC library, "BoB". The results indicated that the individual LTR retrotransposons (copia and gypsy-like) were represented by between 90 and 320 copies in the haploid genome, with only evidence of a single location for the LINE. Sequence analysis of the same elements against genome survey sequence gave estimates of between 60 and 570, but no LINE was found. There was minimal evidence for clustering between any of these retroelements: only half the randomly expected number of BACs hybridized to both LTR-retrotransposon families. Fluorescent in situ hybridization showed that each of the retroelements had a characteristic genomic distribution. Our results suggest there are preferential sites and perhaps control mechanisms for the insertion or excision of different retrotransposon groups.

The diversity of retroelements in diploid and allotetraploid Brassica species.

Using universal PCR primers, some 80 fragments of retroelement reverse transcriptase genes were isolated from 16 accessions of the three diploid and three derived allotetraploid species of Brassica in the triangle of U. Sequence analysis showed that the Ty1/copia and LINE-like elements were distinct, while a third clade could be sub-divided into Ty3/gypsy, Athila and virus-like branches, providing evidence that there are multiple sub-lineages within this group normally considered to be gypsy-like elements in plants. The parsimony trees showed no branches correlating with the known genome relationships for the six diploid and allotetraploid Brassica species, probably because members of the element families were present in the common ancestor of the Brassica and, unlike other repetitive sequences, there is no evidence for genome-wide homogenization, although convergent evolution or horizontal transfer cannot be ruled out. Southern hybridization suggested some sub-families were amplified in individual species. The data show that retroelement sequence data do not allow inference of phylogeny, but knowledge of evolution of such abundant sequences assists in exploitation and interpretation of data from other species including models with much smaller genomes and may provide markers.